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The superconducting critical temperature Tc of a superconductor/ferromagnet (S/F) bilayer with 
spin-flip scatterings at the interface is calculated as a function of the ferromagnet thickness dp in the 
dirty limit employing the Usadel equation. The appropriate boundary conditions from the spin-flip 
scatterings at the S/F interface are derived for the Usadel equation which includes the spin triplet 
pairing components as well as the spin singlet one. The spin-flip processes induce the spin triplet 
pairing components with s-wave in momentum and odd symmetry in frequency from the s-wave 
singlet order parameter A of the superconductor region. The induced triplet components alter the 
singlet order parameter in the superconductor through boundary conditions at the interface and, 
consequently, change the Tc of an S/F bilayer system. The calculated Tcidp), like the case of no 
spin-flips, shows non-monotonic behavior which typically decreases as ciF is increased from and 
shows a shallow minimum and then saturates slowly as dp is further increased. It is well established 
that as the interface resistance (parameterized in terms of 7^) is increased, the Tc is increased for a 
given dp and the non-monotonic feature in Tcidp) is strongly suppressed. As the spin flip scattering 
(parameterized in terms of 7™) is increased, on the other hand, the Tc is also increased for a given 
dp, but the non-monotonic feature in Tc{dp) is less suppressed or even enhanced, through the 
formation of the spin triplet components. 
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I. INTRODUCTION 

The superconductivity and ferromagnetism are two 
competing orders: The former prefers a spin anti-parallel 
state and the latter a spin parallel state. They can coex- 
ist and exhibit interesting interplay effects only in a very 
narrow range of parameters. The effects of their compe- 
tition and interplay, therefore, can be more conveniently 
studied when the interactions responsible for the two or- 
ders are confined to spatially separate regions of a sys- 
tem like a superconductor/ferromagnet (S/F) junction. 
When a superconductor is brought into a contact with 
a ferromagnet, the two competing orders influence each 
other in the vicinity of the interface on a spatial scale of 
the order of the coherence lengths. These phenomena of 
mutual influences are referred to as "proximity effects" . 
The proximity effects of a superconductor/normal metal 
(S/N) junction show up as a penetration of the supercon- 
ducting pairing amplitude into the normal metal region 
with an exponential decay^'^. In an S/F junction, on 
the other hand, the superconducting pairing amplitude 
'i'{x) in the ferromagnet region does not simply decay ex- 
ponentially but also oscillates. This oscillation shows up 
because the electrons forming Cooper pairs feel a distinct 
potential depending on their spin due to the exchange 
field in the ferromagnet, which leads to a finite momen- 
tum of a Cooper pair similar to the "Fulde-Ferrel-Larkin- 
Ovchinikov state" (FFLO state) in bulk materials'^'**. 

The same physics behind this pairing amplitude oscil- 
lation manifests itself in a number of different contexts: 
For instance, the non-monotonic dependence of Tc of 
S/F systems on the ferromagnet thickness dp and the "tt- 



state" of S/F/S junctions, among others. Here, the Tc is 
the superconducting critical temperature of an S/F junc- 
tion where the current is parallel to the interface. Then, 
the Tc corresponds to the highest temperature at which 
the superconducting order parameter A(x) becomes non- 
vanishing at least at one point within the junction. The 
A(x) defined in Eq. (10) and ^{x) are related as 



A{x) ^ A*(a;), 



(1) 



where A is the superconducting pairing interaction. The 
TT-state refers to a case where two superconductors of a 
S/F/S junction separated by a ferromagnet of an appro- 
priate thickness have the phase difference of tt in a ground 
state without any external gauge field. That is, the or- 
der parameters of two superconducting regions have the 
opposite signs. The magnetic coherence length, which de- 
termines the oscillation length and the penetration depth 
of the pairing amplitude ^'(x) in the ferromagnet, the 
ferromagnet thickness dp of minimum Tc in Tc{dp), the 
ferromagnet thickness of the 7r-junction, and so on, is 
given roughly by 



— 



TT-EcX 



for clean limit, 
for dirty limit, 



(2) 



where i?cx is the exchange energy and £ is the mean free 
path. For a ferromagnet, E^^ is about thousands Kelvin, 
which gives a S^')^ of a few nanometers. Studies of the 
oscillatory behavior of S/F junctions, therefore, require 
a fabrication technique of a nanometer scale control of 
the F thickness. 



Due to the difficulties in the fabrication, serious stud- 
ies of S/F junctions were initiated by theoretical works in 
the late 70's. Bulaevskii et al. considered magnetic impu- 
rities in Josephson junction^. Buzdin et al. predicted the 
critical current® and Tq oscillation'' as a function of dp, 
and studied the effects of ferromagnct layer in Josephson 
junction®. Following these theoretical works, many ex- 
perimental efforts have been made to test the pairing am- 
plitude oscillation in the ferromagnct. One of them is to 
measure the Tc vs. dp of S/F bilayers. The measured Tc 
showed a non-monotonic dependence on dp^, although 
there were some contradictory reports'^'^. Another line of 
research is to verify the 7r-statc. Quite a few experiments 
reported the 7r-statc: tunneling spectroscopy-'^^, tempera- 
ture dependence^^, and thickness dependence of the crit- 
ical current^^, and even phase sensitive measurements^^. 
Through all the above experiments, the pairing ampli- 
tude oscillation in the superconductor-ferromagnet het- 
ero structures is now well established. 

The advancement of the experimental technique of the 
S/F systems, on the other hand, demanded more detailed 
theoretical analysis of the systems to include those effects 
hitherto not considered such as the inhomogeneity of the 
magnetization. For example, there were some measure- 
ments that reported much longer penetration lengths in 
S/F jimctions^'"''^®. The Tc difference between the paral- 
lel and anti-parallel magnetizations of the two ferromag- 
nets sandwiching the superconductor in F/S/F junctions 
is smaller by the factor of 10^ than the theory predicts^^. 
For filling the gap between the theories and experiments, 
the spin-orbit scattering in the ferromagnct^*, or the 
triplet components^^"^* began to be considered by sev- 
eral groups. It should be pointed out that there is no 
pairing interaction in the triplet channel in the S or F re- 
gion because conventional s-wave superconductors were 
considered for S/F bilayer junctions. 

The triplet pairing components are induced in a hybrid 
junction of S and F when a single quantization axis for 
spin can not be defined for the whole junction. Such sit- 
uations arise and the triplet pairings were studied for (a) 
the non-coUinear alignment of the magnetic moments of 
two F's in a F/S/F structure22-24 pjg ^ i^^y^^ ^ 

tion of the magnetization direction within a finite thick- 
ness in a F region of S/F bilayers^^ [Fig. 1 (b)], (c) spin- 
orbit scattering near the interface between non-magnetic 
and magnetic materials in S/F structures'*^, and (d) spin 
mixing near the interfaces of S/F/S structure where F is 
a half-metaP^. The works of (a)-(b) considered the dirty 
limit cases. In the dirty limit, the non-s wave compo- 
nents are strongly suppressed, and the so called "odd" 
frequency triplet superconductivity is more robust than 
the p-wave triplet superconductivity. The odd triplet 
pairing has even symmetry in spin and momentum spaces 
but odd symmetry in frequency, and was first proposed 
by Berezinskii'^ in 1974 in the context of '^He superfluid- 
ity. The odd frequency triplet pairings, interestingly, can 
be realized in hybrid junctions of S and F as discussed in 
(a)-(b) which employed the Usadel equation in the dirty 
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FIG. 1: Some schematic pictures of cases which can induce the 
triplet pairing components, (a) non-collinear magnetizations 
of two ferromagnets separated by a superconducting layer, 
(b) rotation of magnetization in a ferromagnct layer, and (c) 
spin-flip scattering at the interface of the magnetic and the 
non-magnetic layers due to mixing of two materials or spin- 
orbit scattering. 



limit'®. The work (c)-(d) considered the p-wave pairing 
in the clean limit. 

In this paper, we considered an S/F bilayer as in (b) 
and (c), and studied the effects of the induced triplet pair- 
ing components in S/F bilayers by employing the Usadel 
equation. The difference is that we consider the spin 
flip scatterings at the interface. Technically, this calls 
for new boundary conditions at the interface of an S/F 
junction. We derived such boundary conditions accom- 
modating the spin-flip scatterings at the interface. The 
work of (b) simplified and bypassed the problem of the 
new boundary conditions by considering an artificial ro- 
tation of the magnetization of a fixed amplitude in the F 
region. On the other hand, we consider an intrinsic mag- 
netic inhomogeneity at the interface of S/F junctions as 
depicted in Fig. 1 (c). The Usadel equation is appropriate 
in the dirty limit and is widely employed for analyzing 
various S/N, S/F, and S/N/F systems since various junc- 
tions of S, N and F films belong to the dirty limit. Here, 
we calculate the critical temperature Tq as a function of 
the ferromagnet thickness dp of S/F bilayers including 
the triplet pairing components in addition to the singlet 
one. 

This paper is organized as follows. In Sec. II, we will 
introduce the notations and write down the Usadel equa- 
tion and boundary conditions that need to be solved. For 
the boundary conditions of Usadel equation with spin-flip 
scatterings, the boundary conditions of the Eilenberger 
equation is flrst derived in the Appendix A and, in Ap- 
pendix B and C, the detailed derivation of the boundary 
conditions for Usadel equation is collected. This will ren- 
der the paper more easily accessible to the less technically 
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inclined. Then, the Usadel equation with the appropri- 
ate boundary conditions is mapped onto an eigenvalue 
problem [See Eq. (24) below.] whose detailed derivation 
is given in Appendix D. Actual calculations of Tc consti- 
tute setting up the matrix Kij of Eq. (24) and finding its 
smallest eigenvalue self-consistently. The self-consistency 
is achieved via iterations. In Sec. Ill, the results of nu- 
merical Tc calculations of S/F bilayers with the spin- flip 
scattering induced triplet components will be presented. 
Overall, the Tc of S/F bilayers with the spin- flip scatter- 
ings, like the more familiar case of no spin-flips, typically 
decreases initially as dp is increased from 0, then shows a 
minimum at a finite dp, then increases slowly and satu- 
rates to Tq = Tcidp oo) as dp is further increased. As 
the spin flip scatterings are increased, compared with the 
corresponding no spin flip case, (a) the Tc is increased 
with decreasing non-monotonic behavior, but (b) the rel- 
ative non-monotonic feature in Tc{dp) with respect to 
T^ is enhanced (the increase of T^ is faster than the 
minimum Tc), through the formation of spin triplet com- 
ponents. Finally, we have summary and some concluding 
remarks in Sec. IV. 



II. FORMULATION OF S/F BILAYER 

In this section, we will present the Usadel equation 
with the induced odd-frequency s-wave triplet pairing 
components in addition to the dominant singlet pairing 
order parameter, boundary condition at the S/F inter- 
face, and calculation of Tc by mapping onto an eigen- 
value problem. The detailed derivations of the boundary 
conditions and mapping of Tc to the eigcnsystem are col- 
lected in the Appendix B, and D, respectively. We con- 
sider the dirty limit case which almost all experiments 
belong to. 

The Usadel equation which considers both the singlet 
and triplet pairing components can be written using the 
anomalous function in spin space as follows: 



^""wTcg^Fix^iuj) = \uj\F(x,iLu) - A{x) 
-I sgn (a;) (jyexF(a;,ia;)-F(a;,ia;)4*,), 



(3) 



where F is the Usadel function from anomalous Green's 
function of 2 x 2 matrix in spin space and A is the singlet 
pairing order parameter. 



F = 












A" 




^11 . 




-A 






(4) 



Hex is also 2x2 matrix which denotes the exchange field 
in the ferromagnet. 

H = H^^a, + Ry^ay + H^^a, (5) 

where Ci's are the Pauli matrices. To show the coupling 
between the singlet and triplet pairing state clearly, the 



following vector form of the Usadel equation can be con- 
venient in numerical calculations. 



Q2 

f'KTc-g^F{x,iw) = \u;\F{x,iw)- A{x) 

-i sgn (w) Hex • F{x,itj), 



(6) 



where lj = 27rT(n -|- 1/2) is the Matsubara frequency, 
D = VpT/3 is the diffusion constant, and 



C = VD/2nTc 



(7) 



is the coherence length in dirty limit. The F and A are 
four component vectors and Hex is a 4 x 4 exchange field 
tensor given, respectively, by 



F{x, iu)) 



H. 



Fs{x,iu!) 
Ftx{x,iuj) 
Ftz{x,iuj) 
Fto{x,iuj) 



A(a;) 



A(x) 








TJX 



^ex -^ex *-^ex 



(8) 



A(.t) is the singlet s-wave pairing order parameter and 
a function of x which represents the coordinate perpen- 
dicular to the interface between S and F. The F^s are 
transformations of the F„„i by expanding F on the basis 
of the Pauli matrices as shown in Eq. (B14) in Appendix 
B. The transformations are given as: 



= -^{F^i- Fi^), Ft^ -- 



:i(i^TT+ni)- (9) 



The Fs is the singlet pairing component and the others 
are the triplet pairing components. The self- consistency 
relation is given by 



/^{x)=\-KTcg{ep)Y,F, 



(10) 



where A is the superconducting pairing interaction and 
giep) is the density of states per spin at the Fermi level. 
Since we assume that A is zero at ferromagnet, we should 
note that the superconducting order A is also zero at 
ferromagnet and that only the pairing amplitude $ can 
penetrate into. 

The Usadel equation of Eq. (3) or (6) must be supple- 
mented by a set of appropriate boundary conditions. The 
system we consider is an S/F junction where the inter- 
face between S and F has the potential (without spin-flip) 
and spin-flip scatterings due to magnetic inhomogeneity 
as shown in Fig. 2. Therefore, the interface potential is 
represented by 6 functions with the spin dependency as 

U{x) = {Voao + V^a^ + VyCFy + V^a^) 6{x). (11) 
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For simplicity, we can set Vq and T4 to be the only non- 
zero parameters among the spin dependent 5-function, 
when the magnetization of ferromagnetic layer is in z- 
direction. The spin-flip scatterings were modeled in 
terms of the (5-function- like magnetization at the inter- 
face. In real situations, this will be local inhomogeneities 
of the magnetization within a very short range near the 
interface. Within the Usadel formulation, the boundary 
conditions using the similar notation with Eq. (6) can be 
written as 

F^-F^ = [7b-z sgnH7™]Cj.^F^, (12) 



ox ox 



0, (13) 

where the dimensionless parameters are defined by 



7 



Ps^s 

PF^F 



lb = 



7m = 



76 + 71 

76 - 71 



72 
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7m 








76 + 71 
76 





72 





71 



(14) 



The 7 is the ratio of S and F bulk properties, and 7;, rep- 
resents the contact quality which is expressed as 
without spin-flip scatterings. A and TZ are the inter- 
face area and resistance, respectively. 7™ is the spin-flip 
strengths between S and F which is proportional to Vx 
in (11). 71 and 72 comes from the higher order term 
of Vx and Vq which can be ignored when the scattering 
strength is weak. See Appendix B for details. 
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FIG. 2: The hypothetical layer (I) with an infinitesimal thick- 
ness 25 between superconductor (S) and ferromagnet (F) 
which are described as a spin-dependent potential U (x) which 
will give a interface resistance and a spin-flip scattering, i-i(x) 
stands for a local chemical potential. 

The Usadel equation of Eq. (6) with the boundary 
conditions of Eqs. (12) and (13) was mapped onto an 
eigenvalue problem and solved numerically to determine 
the Tc- The Tc defined here can be determined exper- 
imentally by measuring the resistance with the current 
parallel to the interface. For simplicity, we take the mag- 
netization at the interface along the x-axis and that of 
ferromagnet along the z-axis. 



We follow Fominov et a/.^^'^^ to get the numerical solu- 
tion of the Usadel equation. We first separate the anoma- 
lous Green's function Fa{x,iLL!) into the even and odd 
symmetry parts in the frequency space 



F^'^\x,iuj) = Fa{x,iu!) ± Fa{x,-iuj), 



(15) 



where a = s, tx, tz, and tO. Since all the anomalous 
Usadel functions are s-wave pairing, we choose the even 
symmetry in frequency for the singlet component, and 
the odd symmetry for the triplet components. Thus we 
define a 4-component vector as 



^(x, iuj) 



F^+\x,iuj) 



tz 

A') 



(x, iuj) 
(x, iuj) 
(x, iuj) 



(16) 



Here, the triplet components, Ftx, Ftz and Fto, are odd 
in the frequency and even in the momentum space as we 
mentioned in Sec. I. This special type of triplet pairing 
is the odd triplet pairing proposed by Berezinskii^^. It 
will be much more robust than a p-wave triplet pairing to 
disorders due to the momentum independent gap feature. 

As shown in Fig. 3, for an S/F bilayer, we took the 
superconductor region in < a; < dg, the ferromagnet in 
—dp < X < 0, and the spin-flip scatterings at the inter- 
face oi X — 0. In the S region, there is no exchange field 
and only pairing order parameter A(a;) exists. Therefore, 
the Usadel equation of (6) can be written for the S region 
as 

TrTces^^H^, = H^'^i^, i^) - 2A(x), (17) 

for < X < ds- For the F region, on the other hand, 
there is no pairing interaction and only the exchange field 
term appears. The Usadel equation is then written as 



'^^cip-Q^^^ {^1^^) = (kl - iHox)^^(a;,iw), (18) 



for —dp < X < 0. The coherence length ^p in the F 
region is given by Eq. (7) as 
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2ttTc 



(19) 



where Dp is the diffusion constant of F and Tc is the 
superconducting transition temperature of the S/F junc- 
tion. The £^p should be distinguished from the ^™ of Eq. 
(2). S^f!^ corresponds to the actual penetration depth of 
singlet pairing amplitude in F, but is the pure pene- 
tration depth without considering the exchange field as 
in normal metal. Later, this ^p will turn out to be a pen- 
etration depth of triplet components, Ftz and Fto, which 
is unaffected by the exchange energy. The boundary con- 
ditions of Eqs. (12) and (13) can be written in terms of 
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J-{x,iuj) of Eq. (16) as follows: 



III. NUMERICAL RESULTS 
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(20) 
(21) 
(22) 

(23) 



Here, 71 and 72 of Eq. (14) are ignored for the weak 
scattermg rate at the interface. 




FIG. 3: A schematic plot of the pairing amplitude $(x) of 
Eq. (10) with the dotted line for an S/F bilayer. The magne- 
tization in the ferromagnet region is taken in the z-direction 
and the magnetization at the S and F interface is in the x- 
direction. 

As detailed in Appendix D, Eqs. (17), (18), (20), (21), 
and (22) can be solved for the Tc of a S/F bilayer, and 
the Tc equation with respect to the bulk transition tem- 
perature Tco may be cast into the form 



Ailn 



Tco 



N 



(24) 



where the order parameter A and kernel K are defined 
on a discrete grid as 



A, = A I'^z 

N 



2vrTc J2 
ds 



- G.ij{iuJn) 



(25) 



Here, N is the number of the discrete grids in the inte- 
gral domain between and ds- The Tc is determined 
by the smallest eigenvalue of the matrix K given by 
(25), which gives the largest Tc- However, the kernel 
Kij is Tc dependent through the Matsubara frequencies 
ujn = 2TTTc{n + 1/2). Therefore, we solve the Eq. (24) 
iteratively until it gives a self-consistent solution for Tc- 



The interface parameter of the boundary condition, 7f, 
in (21), contains diagonal and off-diagonal elements as 
given in (23). The diagonal element, 7;,, corresponds to 
the interface resistance and the off-diagonal element, 7^, 
to the scattering between singlet and triplet pairing com- 
ponents. The relation between the parameters 7f,, 7„i and 
the more microscopic processes like the potentials Vb, Vx 
is complicated and depends on detailed mechanism. In 
the present work, therefore, we take the 7^ and 7,„ in 
(23) as independent parameters characterizing the S/F 
interface without explicit references to more microscopic 
processes. For the present microscopic consideration in 
terms of Vb and Vx , the connection between the interface 
parameters and the microscopic potentials can be found 
in Appendix B. 

The Usadel equations in Eqs. (17) and (18) with the 
boundary condition in Eqs. (20), (21), and (22), can be 
expressed by the seven dimensionless parameters, ds/^s, 
dp/^F, Tc/Tco, Hc^/Tco, 1, lb, and 7„. All the calcula- 
tions for Tc /Tco are done and shown in the figures with 
respect to the other six independent parameters. The nu- 
merical results were calculated based on the parameters 
in Ref. 28, which are Tco = 7.0 K, is = 8.9 nm, ^f = 7.6 
nm, ps = 7.5 /ifJcm, pp = 60.0 /ificm, and ds = 11 nm. 




FIG. 4: Tc /Tco vs. cIf/^f for the various Hex without in- 
terface resistance and spin-flip scattering at the interface, i.e. 
76 = 0.0 and 7,„ = 0.0. 

First, we show the typical behavior of Tcidp) in S/F 
bilayer in Fig. 4 as well as Tcid^) in S/N bilayer 
(solid line) which is S/F bilayer with zero exchange field. 
Hex = 0.0. For S/N bilayer, Tc decreases monotoni- 
cally as the thickness of normal metal layer increases. 
When the exchange field is turned on, the Tc's, the dot- 
ted curves in Fig. 4, show typical behaviors of Tcidp) 
for S/F bilayers. As dp is increased from 0, Tc decreases 
fast initially and shows a minimum, and then increases 
slightly and saturates as dp is further increased. This 
non-monotonic dependence of Tc on dp is due to the os- 
cillation of the pairing amplitude in ferromagnet, known 
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as FFLO state. The thickness of minimum Tc in Tc{dp) 
is almost corresponding to the dirty hmit value of in 
Eq. (2). For the same reason, the minimum Tc can be 
also denoted as the transition from 0-state to 7r-state of 
pairing amplitude in F. As shown in Eq. (2), the posi- 
tion of minimum Tc is shifted to left as If ex increases. 
For a strong enough exchange energy, Tc is strongly sup- 
pressed and goes to zero as dp increases. For a properly 
strong exchange field, moreover, the reentrant behavior 
of superconductivity occurs as the cases of ifex = 120 K 
and 150 K shown in Fig. 4. 



1.2 

1 

0.8 



,0.6 
0.4 
0.2 




7^=0.0 

7^=0.05 

\,„ Y6=1.0 













0.2 0.4 0.6 0.8 



1 1.2 1.4 



In Fig. 6, Tcidp) is calculated for various 7™, the 
interface spin flip parameter, with the other parameters 
same as in Fig. 5. A non-zero 7™ induces the triplet 
pairing components and couples them with the singlet 
component. As can be seen from the figure, Tc is en- 
hanced by jm as well. Recall that Tc is determined by 
the singlet pairing order parameter A (a;) in the S region, 
and there is no pairing interaction in the triplet chan- 
nel for the S/F system we consider here. The triplet 
components come into existence only through 7^, and 
an increase of 7^ means an increase of the triplet pair- 
ing components. The Tc is enhanced because the singlet 
paring order parameter A (a;) in the S region is enhanced 
when it is coupled to the triplet components in addition 
to the singlet pairing component. The exchange field of F 
region does not act as a pair breaker for the triplet com- 
ponents unlike for the singlet one. The detail effects of 
the 7to is somewhat different from 7;,. When we compare 
two cases of 75 = 0.2 in Fig. 5 and 7™ = 0.5 in Fig. 6, 
both have almost same saturated T^ = lim^p^oo Tc{dp), 
but the non-monotonic structure is much stronger in the 
later case than the former one. While the 7^ is similar to 
7fc in the sense of the increase of Tc and the suppression 
of a minimum Tc at a particular dp, they are different 
in that the non-monotonic structure is less suppressed 
by 7to than by 75. This can be more clearly when the 



FIG. 5: Tc/Tco vs. dp/^p for the various 75 without spin-flip 
scattering at the interface, i.e. 7^ = 0.0. 
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FIG. 6: Tc/TcQ vs. dp/^p for the various 7„ 
scattering at the interface with "ft = 0.0. 
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FIG. 7: {Tc-Tc)/{Tco-Tc) vs. dp/^p for several values of 
the interface resistance parameter "ft with the interface spin- 
flip scattering parameter 7^, = 0.2. As the jb is increased, 
the non-monotonic feature is successively weakened. The Tc 

increases as 75 is increased but appears the other way because 
it is normalized as Eq. (26). 



Fig. 5 shows the effects of interface resistance on the 
F thickness, dp, dependence of the superconducting crit- 
ical temperature Tc of S/F bilayers. All the parameters 
except for 7;, are same as used in Fig. 4 and the exchange 
field is fixed to i?ex = 150 K. As the interface resistance 
parameter 7;, is increased, the Tc{dp) is increased and 
the non-monotonic feature becomes weakened. The Tc 
is increased because a large 76 decouples S from F and, 
therefore, S is not influenced by F. Fig. 5 shows that the 
non-monotonic feature is almost washed out for 7^ > 0.2. 



Tcidp) is normalized in the following way: 

{Tc-T5)/{Tco-m. (26) 

Fig. 7 and Fig. 8 show how the non-monotonic fea- 
ture changes as 7^ and 7^, respectively, arc increased 
from 7b = 0.05 and jm = 0.2. Both figures were calcu- 
lated with the same parameter set as in Fig. 5 and Fig. 
6 except for the 7^ and 7^. The results for = 0.05 
and 7to = 0.2 are common for both figures and denoted 
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FIG. 8: {Tc-Tc)/{Tco-T^) vs. dp/iF for various 7^ with 
76 = 0.05. Notice that an increase of 7™, enhances the non- 
monotonic feature. 



as the solid lines. In Fig. 7, the minimum values of the 

normalized Tc increase, or the non-monotonic feature 
is weakened, as 7t increases, vifhile the minimum values 
decrease due to 7™ in Fig. 8. This means that the en- 
hancement of minimum Tc is less strong than that of Tq. 
The reason is that the magnitude of the induced triplet 
components depend on the thickness of the ferromagnet 
of a S/F bilayer. The length scale over which the triplet 
components develop is given by ~ l/k^p as in (D6) 
which is longer than that of singlet component in ferro- 
magnet given by of Eq. (2). Therefore, the Tcidp) 
with non-zero jm, compared with Tcidp) with 7^ = 0, 
is increased more for dp > £,f than for dp < £,p. 

Another way of seeing why the magnitude of the in- 
duced triplet components depend on the ferromagnet 
thickness is from the boundary conditions of Eqs. (20) 
and (21). The boundary condition of (21) means that 
the magnitude of the triplet components are proportional 
to the derivative of the singlet component at the inter- 
face. Since the derivative of the pairing components at 
X =- —dp is equal to from (20), the derivative at .x = 
will be close to when dp approaches 0. This observation 
means that the magnitude of the triplet components also 
go to as dp goes to 0, and the enhancement of Tc due to 
the induced triplet components is not as strong for small 
dp &•& for large dp. Since the di? at which the Tc{dp) is 
the minimum is about which is smaller than ^p, the 
enhancement of Tc for dp > ^p gives relative increase of 
the non-monotonic structure. This enhancement of the 
normalized non-monotonic feature is quite distinguish- 
able from the effects of interface resistance and spin-orbit 
scattering in bulk ferromagnet^*, which give the suppres- 
sion of the non-monotonic structure. 



IV. CONCLUSION 

In this paper, we considered the effects of the odd- 
frequency s-wave triplet pairing components on the tran- 



sition temperature vs. the ferromagnet thickness of S/F 
bilayers. The triplet pairing components are induced by 
the spin-flip scatterings at the interface. The interface 
spin-flips physically occur because the magnetic and non- 
magnetic atoms diffuse around the interface between the 
magnetic and non-magnetic layers, because the direction 
of magnetization is changed locally near the interface due 
to the boundary effects, or because the spin-orbit scat- 
tering at the interface induced by the polarization due 
to the difference of the work-function of two adjacent 
material. The spin-flip scatterings at the interface were 
simply modeled in terms of a spin-dependent (5-function 
which includes potential scattering and spin-flip scatter- 
ing. The appropriate boundary conditions describing the 
spin-flip scatterings at the interface between S and F were 
derived in the context of the Usadel equation and Eilen- 
berger equation. Through the boundary condition, the 
singlet and triplet pairing components are coupled each 
other. The Usadel equation with the derived boundary 
conditions for S/F bilayers was then mapped onto an 
eigenvalue problem. It was solved self-consistently via 
iterations to obtain the superconducting transition tem- 
perature Tc of S/F bilayers. 

Compared with no spin-flip cases, Tc is enhanced by 
inducing the triplet pairing components via spin-flip scat- 
terings at the interface. In detail, this enhancement is 
distinguishable from the enhancement of other effects, 
for example, the interface resistance and the spin-orbit 
scattering in bulk ferromagnet. All these effects mostly 
suppress non-monotonic dependence in Tc vs. dp, but 
the enhancement of Tc due to the spin flip scatterings is 
accompanied by an enhancement of the relative depth 
of minimum Tc which is scaled by the difference be- 
tween Tco and the saturated value, Tq. This enhance- 
ment of the relative non-monotonic structure in Tc{dp) 
is related to the coherence length S^p which is the decay 
length of triplet pairing state in ferromagnet. But the 
triplet paring components affect Tc{dp), and can be dis- 
tinguished from other effects by quantitative analysis of 
Tc{dp) measurements or tunneling experiments. 

Wc arc currently extending the present analysis to the 
superconductor/normal metal/ferromagnet (S/N/F) tri- 
layer systems. The results will be reported separately. 
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APPENDIX A: BOUNDARY CONDITIONS FOR 
EILENBERGER EQUATION 

Since the quasi-classical equations such as Usadel equa- 
tion and Eilenbcrgcr equation describe variation of phys- 
ical quantity over longer length scale than interatomic 
length scale, they are invalid in dealing with the effects 
of interface directly. To include the spin-flip effects at the 
interface in the quasi-classical equations, we will impose 
such effects on the boundary condition. Since the Usadel 
equation is originated from the Eilenberger equation and 
the Eilenberger equation is again from Gor'kov's equa- 
tion, the derivation of the boundary condition consists 
of two steps. First, we will derive the boundary condi- 
tion for Eilenberger equation which is the quasi-classical 



approximation of Gor'kov's equation. From Eilenberger 
equation, we will derive the boundary condition for Us- 
adel equation which is the dirty limit of Eilenberger equa- 
tion. The boundary condition for the Eilenberger equa- 
tion without any spin-flip scattering at the interface is 
originally derived by Zaitsev^^. Therefore, we will fol- 
low his procedures and then include the spin degrees of 
freedom and spin-dependent interface potential for the 
spin-flip scattering at the interface. For the insufficient 
part in this derivation, readers refer to Ref. 29. 

First, we write down the Gor'kov's equation, which 
describes the system of two adjacent region separated by 
intermediate region with thickness 25, as shown in Fig. 
2. 



(5(r-r')5(T-T')i= (-p.|: + ^^ + A-C/(r)-E + A^(r)) G(r, r'; r, r'), (Al) 

I 



where G is 4 x 4 Green's function in Nambu-Gor'kov 
space. /z(r) is the chemical potential, and E is the self- 
energy. Here the U{y) is the interface potential including 
spin-flip scattering, which is expressed as the following 
spin dependent delta function in spin and particle-hole 
space. 

C/(r) = 5{z) (VoPoo-Q + VxPoOx + Vyp^Oy + V^poaz) . 

(A2) 

Pi and (7i is Pauli matrix in particle-hole space, and in 
spin space, respectively. For the convenience, we consider 
Vo and only for potential scattering and spin-flip scat- 
tering, respectively, and ignore the other terms without 
loss of generality. Now we separate the above Green's 



function into fast varying part and slow varying part on 
z. 

G{z,z';p,ioJn) = Gn{z,z';p,iujn)e'P^'^'-'"> (A3) 

+Gi2(z,z';p,*w„)e^f^(^+^') 
+G2i{z,z';p,iiVn)e-'P'^'+''^ 

Now Gij is the slowly varying function of z and z' . If we 
substitute Eq. (Al) by Eq. (A3) and neglect the second 
order derivatives, we get the equation for Gij . 



d 



d 



iPzOJn - (-1) Wzi ^ + ttVii ^ + A(r) - E 



dz 



dp 



Gkn{z, z'; p, iUn) =0, {z^ z') 



(A4) 



r 



Here, we assume the z-direction is normal to the interface which are continuous at z = z'. 
of the two region and p is the coordinates normal to the 
z-direction. We assume the translational symmetry in 
p-direction, ans also assume that the boundary is sharp (l(z,z ) = 2i\v 
and plane. When z,z' > 0, v^i is replaced by Vz2- The 
function Gij 's for z < z' and z > z' are matched by 



Q{z,z') 



2iu 



Gii{z,z') - sgn (z - z'), p^j > 
G22(z,z') + sgn (z - z'), Pzj < 

Gi2(z,z'), > 0, 

G2l(z,z'), Pzj < 0. 



(A6) 



Gfc„(z'+0,z')-Gfe„(z'-0,z') = (-1)'=- 



-Skn- (A5) 



To make the equations simpler, we introduce g and G 



For z <C z* = min(^j, vpj /e, a), we can write the equation 
for the Green's function using the fact that ^ ^ ^ ^ 
eG, where e ~ max(T, A, V, ui) as done by Zaitsev^^. The 
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equation for Green's function G is reduced to 
H{z)G{z, z') = H{z')G{z, z') = 6, 

where 



Pii 



(A7) 



(A8) 



for z ^ z'. The solution can be written as 

G = ^^e'^'^^"^'' +^±^e-'f'(^-^') (A9) 

+B(j)e'P«(^+^') + BO)e'P-(^+^')^ (^zz' > 0) 
G = A(-'''='e'*'^(^-^'^ + A(-'''='e-'P»(^-^'^ (AlO) 

where the superscript j, k = 1(2) denote the index of the 
region z < —5 {z > S), and the index +(— ) eorresponds 
to z > z' {z < z'). Now the above Green's functions can 
be written down by using the two hnearly independent 
solution of the Schrodinger equation H{z)il)i^2{z) = 0. 



V'2 



ptp.iz j^fQ-iPzlz^ (Z < 0) 

de'P-^'', (z > 0) 



(All) 



H^de^P-^^ (z < 0) 

- ^ ' (A12) 



ft^e-*P-l^ {z > 0) 



where f and d is the spin dependent coefficient. For the 
spin space, it will have the form, 



(Pzi+P.2+iVo)^ + K^' (p.i+P.2+iVo)^ + V^ 
-2ip,iV^ pli-V^-{p.2+m) 



-2ip^iV^ 



2pzl(Pzl+P»2+»V"o) 

(p.i+Pz2+iyo)2 + K„2 (p,i+p,2+iVh)2 + V2 

-2ip^iVa; 2p^i(p^i+p^2+»Vb) 



,(A13) 
.(A14) 



For 1^1, \z'\ <^ z* , G can be represented in the form, 

G{z,z')= (A15) 
' e^P-^^'Fi{z) + e-'P-i^'Fi(^), -2' < z> z' 



e»p.2^^2(z) + e- 



tj»P2l2 



Pi(z') 



-JPzlZ 



Pi(z'), z < -5, z < 2' 



where 

P,(z) = V'i(^)-fa+V'2(^)-fi2, (A16) 
Pi (2;) = il)x{z)-Vii^i>i{z)-Vii, (A17) 

Here, fii(2) and Pii(2) are 4x4 matrix in spin and particle- 
hole space. The coefficients , B^^ and B^^ is 
related to the Green's fimction at the vicinity of the in- 
terface. For the relation between the Green's function of 
the region 1 and region 2 near the interface, we compare 
the solution of Eq. (A15) with Eqs. (A9) and (AlO), then 
we can write down the coefficients B\ B^^ and B^^ 



in terms of fii(2) and Pii(2). If we eliminate the term 
fji(2) and Pji(2), we get the following relations. 



vxAS^ - V2A% = viA]^ - V2A% = 

- r+V2^' 
d _ _dt 

UlA^-hV2^± = ViA]^+V2A% = 



V2^B'^ +vif^B^ = vifB^ +V2'-^B'^, (A18) 



-y^^^B^ + y^l^B' =-V2^B^ + v,Ib\A19) 
Here, we define the function Q from Q. 



rQ z < —S, Pi > 

r^g z < -5, Pi <0 

S=< r^g z>S, pi>0 

i^g z>S,pi<0 



(A20) 



Now we define the symmetric and asymmetric component 
of the Green's function. 



9s(a) 



9{P) ± 9{-p) 



(A21) 



From the Eqs. (A9) and (AlO), we get the following 
expression for g^(^a) and ^^(a): 

g, = ivj (^i + 

9a = iVj - 1±) - sgn (z - z') 

( ^{fB^+r^B^) z,z'<-5 

e« = I ^ j^r|t^2 + ^'52^ Z,Z'>6 



ga = 



i^{fB^~r^B^) z,z'<-6 



(A22) 



From the relation Eqs. (A18), (A19), and (A22), we 
get the relations for the symmetric and the asymmetric 
Green's function. 



9i'^ = Si^) = 9a 

gi'^ - gf^ - ga 



(A23) 
(A24) 
(A25) 



e^+e^ = ^(gi'^-gi'^) (A26) 

Here, R is refiection coefficient. The refiection and trans- 
mission coefficient are 

D = ^dU, R = fH = l-D. (A27) 

Pzl 

For the derivation of the boundary conditions in terms 
of g only, we use the following relation. 



gg = {-lysg-apj g 



(A28) 
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For the proof of the above relation, the readers refer to 
Zaitsev^^. By using (A28), we get the foUowing relations. 



9sGs +9aQa = {-lyQa, 
9eQa+9aQs = {-lYQs- 



(A29) 



Using the above equations with Eqs. (A23), (A24), 
(A25), and (A26), we eliminate Qs and get the follow- 
ing relations. 



-9tR9t - {9s y 

9t9s +9sR9t 

gtSa - gags 
gjQa - gaRgt 



Ga, 
gaGa, 

-Rgt, 
-97, 



where we introduce 



(A30) 
(A31) 
(A32) 
(A33) 



(A34) 



From these four relations, we again eliminate Ga and we 
get 

gt {g^f + ga {g7)'= {i - (gt)'} Rgfg^ (A35) 
[g^f gt + gaR {gtf = g^gt { i - ^ (5+) ' }(A36) 

By considering the following commutation relations, 

[5^,57]+ =0, and [gf,ga]+=Q, (A37) 

finally, with Eq. (A23), we get the following two bound- 
ary conditions in terms of g only. 



gf' gi'' 



ga 



R{gt)' + {g;y 



0, (A38) 

Dg;gt (A39) 
,2" 



Rgs gt, {gty 



This is the new boundary condition for the Eilenberger 
equation with spin-flip scattering at the interface. 



APPENDIX B: BOUNDARY CONDITIONS FOR 
USADEL EQUATION 

For the case without spin-flip scattering at the inter- 
face, The boundary conditions for the Usadel equation 
is derived from the Eilenberger equations by Kuprianov 
et al.^^ in the vicinity of interface within the mean free 
path £. We will use the same methods and assumptions of 
Kuprianov et al. for deriving the boundary condition of 
Usadel equation from that of Eilenberger equation which 
we derived in Appendix A for spin-flip scattering at the 
interface . The boundary condition is directly connected 
to the continuity of the normal and anomalous Green's 



functions of Eilenberger equation. In the Eilenberger 
equation, if we neglect terms of lo and A, we get 

2£,-Vgi = gl{gl)-{gt)g: (Bl) 
2£,-\7gl = gl{gD-{gDg: (B2) 

where g is the quasi-classical Green's function in Nambu- 
Gor'kov space. 



27riV(0) 



■ig f 
P -ig 



(B3) 



with normal and anomalous Green's function in spin 
space. 



5TT gn 
5iT gii 



, and / = 



/tt hi 
fu fll 



(B4) 



Here, the brackets imply an integration over the solid 
angles, (...) — J dVl/A-K. The superscript i in Eqs. (Bl) 
and (B2) stands for the region index as j and k in Eqs. 
(A9) and (AlO). Far from the boundary, the isotropic 
Usadel function^® is identified as g^, and ga can be also 
written in terms of Usadel function. 



gl = {g\) = Gu gl = £-{GiVGi 



where 



G 



iO F 

Ft -id 



(B5) 



(B6) 



The derivation of the boundary condition of the Us- 
adel equation for the derivatives of the Usadel function 
from the first boundary condition of Eilenberger equa- 
tion, (A38), is the exactly same as that of Kuprianov et 
al.^'^ . Therefore, we get the following boundary condition 
for the derivative of the anomalous function in Usadel 
equation. 

For the second boundary condition, we assume i? - 1 ^ 
D{ggy, and R{ga)^ is neglected. From Eq. (B5), we have 
ffa ^ gl enough from the boundary. With the help of 
the relation -|- = 1, we simplify the Eq. (A39). 



gaR = gs gfD 



(B8) 



From the relation between Usadel function and Eilen- 
berger function, Eq. (B5), we rewrite Eq. (B8) in terms 
of Usadel function. 

4G2VG2^^^^ = G2Gi-GiG2 (B9) 

For the boundary condition of anomalous Usadel func- 
tion, we write the Eq. (B9) in terms of normal and 
anomalous Usadel function by using Eq. (B6). 



[G2, Fi — F2]+ + [Gi, Fi — ^2] + 

(BIO) 
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At Tcj. the normal Usadcl function can be expressed as 
the fohowing form which is derived in Appendix C. 

Gi(2) = (i + Re (f)) sgn {w) + i\m (f) , (Bll) 

Now we substitute (Bll) into (BIO), the boimdary con- 
dition for the anomalous Green's function will be 

^2(1 + »^ + isgn (a;)Im f)VF2^B = 

F1-F2 + [Re f + isgn wim f, Fi - ^^2]+, (B12) 



where 



_ 2^2 /X2R\ 

^^~36 \ D I 



(B13) 



Here the anomalous Usadel function is 2 x 2 matrix in 
the spin space. This can be represented in the form of 
spin-singlet and triplet components. 

F(a;) = ioyF,{-V) + o^Ft^{-) + a,Fu{-) + CToFto(-) 

(B14) 

where 



Ti 



IT 



If wc compare the result of Kuprianov when r^; = 0, we 

get 



70 = 7 = 76 



(B22) 



When the interface potential Vq and are small, and if 
we ignore the higher terms of Vq and Vx , the parameters 
in terms of the lowest order of Vq and V^. 

72 ~« X ='^"^' 74«0(V;2) (B23) 

Pfi + PF2 

Finally, the two boundary condition for the Usadel func- 
tions are 



6V:Fi = 76V:F2, 7 = 



IT 



P26 

7fc i^rn 

76 

«7m 76 

. 76 



(B24) 



- -f^TT 



ii 



and Fto = 



i^tt -|- F\ 



- -^TT 



ii 



(B15) 



and the plus and minus sign means even and odd function 
in the frequency space. 

p^^^^ ^ F^{u;)±F„{-uj) ^ ^ 

(B16) 

In vector notation of the anomalous Usadel function we 
have the following form of the second boundary condi- 
tion. 



where 



FU-) 
Fto(-) 



, and 76 



71 72 
73 74 

72 7i 
74 73 



(B17) 



(B18) 



Since we expressed the spin-dependent interfacial poten- 
tial in terms of ctq and a^, the following quantities are 
also expressed as 



/2X2R\ 

(r) = roao + r^a^, ( 3 ) = ^octo + 7xC^x 
Therefore the 7i's in (B18) can be expressed as 

1 / 7a:Rera: \ 70 

^^ = 2r"TTR^J' = 

i 7a;Imro - 7oImr 

X ^x 

^^~2 1 + Rero ' 



(B19) 

(B20) 

(B21) 



APPENDIX C: DERIVATION FOR THE 
NORMAL USADEL FUNCTION AT Tc 

In this section, we will derive the Usadel function of 
normal Green's function at Tc G of Eq. (Bll) in Ap- 
pendix B. Without any spin- flip scattering at the inter- 
face, G is simply as following. 



G = sgn uj for Vx 



0. 



(CI) 



To get the expression for G with spin-flip scattering, Vx ^ 
0. We start from the eigen function of the Hamiltonian 

(A8), 

1 1 r (e'^- + fe-^^-) (x<0) 

^~V2^ IrfV^^^ (x>0) ' (^^^ 



cl>_{x) = 



\JV\IV2 I pi"^ 



{x < 0) 



I (^e-^P^^-ftJ-e'P^^) (a;>0) 



(C3) 



Now, we derive the normal Usadel function from the fol- 
lowing definition of the Green's function. 



G 



N{0) 



E 



iuj — e 



(e>0) 
(e<0) 



(C4) 

Here A^(0) is the density of state at Fermi level. For the 
simplicity, we will deal with one-dimensional case. If we 
write down in terms of the center of mass and relative 
coordinate, we have the following form: 
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i ^ ei{P+p)z + -i{P+p)z I AUi{P+p)Z I ^g-i(P+p)Z 

o = ^i:£ if ±Li ±12 ,(z<.<-zw-z)m 

-i(P+p)z , „i(P+p)z I -U^(P+p)Z , i(P+p)Z 

+^ — ^^^^ e{z <z< -z)e{-z)e{-e) 

f.-i{P-p)z ^ gi{P-p)z _ J.gi(P-p)X _ filg-i(P-p)Z 

+— ^ ei-z <z< z)e(z)e(e) 

P2 IW — € 

(AiP-p)^ + g-i{P-p)z _ f.\^^i(P-p)Z _ filg-iiP-p)Z 

+— ^ o(-z <z< z)e(z)e(-e) 

P2 ILO — e 

3 -i{pZ+Pz) I 3Ui(pZ+Pz) 

+— . {0{Z)e{z < -Z) + 0{-Z)9{z < Z)) 

lu) — e 

3 -i{pZ-Pz) , 3\ i(pZ-Pz) 

+— {e{z)e{z > z) + e{-z)e{z > -z)) 

luj — e 



At the boundary Z = 0, we get 

n"! = , / ..V. ^ , ..... ^ 



j /"OC iPz _|_ „-iPz ^piP^ — p-iPz 

G{k^,uj,Z = = dzY^Red ^-—^ +i\md — — sgn ^ (C5) 

J —oo r> 



After integration over z and P, we have the normal Us- 
adel function with spin-flip scattering at the interface, 

G{k^,ij, Z = 0) = (1 + Re (f))sgn (w) + ilm (f), (C6) 

Without any reflection at the interface, Eq. (C6) is re- 
duced to to Eq. (CI). The imaginary part in Eq. (C6) 
has key role to give a coupling term between singlet pair- 
ing state of even function in frequency and triplet pairing 
state of odd function in frequency. 



APPENDIX D: CALCULATION OF Tc 

In this section, we will get the Tc equation of a S/F 
bilayer by solving Eqs. (17), (18), (20), (21), and (22). 
The triplet pairing components of J-^ of Eq. (18) sat- 
isfy the homogeneous equation unlike the singlet pairing 
component which has inhomogeneous term A (a;). There- 
fore, the solutions of triplet components in Eq. (17) with 
the boundary condition (20) can be written down easily 
as follows. 

Ff^{x) = Cf^ cosh.ks{x - ds), {a = x, z, and 0) (Dl) 
with 



ks 




V ""^c 



(D2) 



In the ferromagnetic region, Eq. (18) is also homogeneous 
equation, and easily solved. The solutions of in Eq. 



(18) with the boundary condition (20) are 



J^^ix) = 



where 



with 



cpix) 
-cf{x) 










4(x) 























'^tz 




r<F 



cf{x) = coshkrix + dp), 
c%{x) = coshkp{x + dp), 



kp = 




andfcp 




(D3) 

(D4) 
(D5) 



(D6) 



Ff , the singlet pairing component of the vector , 
satisfies the inhomogeneous equation (17), but the oth- 
ers, the triplet pairing components, do the homogeneous 
equation. So, it is useful to define projection operators 
to distinguish the two group. 



(D7) 



Here, the subscript s denotes as singlet, and t as triplet. 
At x = 0, from the solutions (Dl) and (D3), we get the 
following relations. 



- 1 








0" 




- 








0" 














, Vt = 





1 


























1 





.0 








0. 




.0 








1 . 



Vt^s^:F^{0) = -AsVtJ^^iO), and 



' dx 



(D8) 
(D9) 
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respectively, with 

As = fcsCs tanh fc^rfs, and 

Re^F -ilmBp 
-ilmBp ReBp 
Bpo 

Bfo 



where 



Bf 



(DIO) 

(DU) 



Bp = [kp^F ta.nh. hrdp] ^, and 
Bpo = [k'^p^F tanh k^dp]'^ . 



where 

(D12) 
(D13) 

Prom Eqs. (D9), (21), and (22), we get the coupled 
boundary condition for J^^ at a; = 0. 

d_ 

' dx 

7/ 



where 



A = 



-Ih 



(D14) 



(D15) 



Since Vs + Vt is identity operator, the following should 
be satisfied. 

{Vs + ^t)?s^^'^(0) = {Vs + Vt)A{Vs + Pt)^'^(O)- 

(D16) 

By utilizing (D8) and (D14), wc can derive the bound- 
ary relation between ^-Ff (0) and F^{0). The boundary 
condition for the singlet component at the S region is 



where 
W 



^s^Ff (0) = W{u;n)Ff{0), 



7[7 + ^s(Re-Bjr+7i)] 



(D17) 



As\Bp + 76|2 + 7(ReBF + 76) + ^S7m«' 
7 + As(ReBF + 7b) 



. ^ • (D18) 

Here, we reduce the three component inhomogeneous dif- 
ference equation to the single component inhomogeneous 
differential equation for F/(x, iiVn), which satisfies (D17) 
and 



g2 



-2A{x) 

(D19) 
from Eq. (17). 

To solve the inhomogeneous equation, we have to solve 
the following source equation. 

^Tce^G{x, y) - \oj\G{x, y) = -5{x - y). (D20) 
ox^ 

With the boundary condition (D17) and (20), we have 



^s^G{Q,y) = W(a;„)G(0,2/), 
d 



(D21) 



0. (D22) 

The solution oFEqs. (D20), (D21), and (D22) has the 
following form. 

ks/\'^n\ 



G{x,y,iwn) = 



smhksds 



ksis 



cosh ksds 



v,{x)v2{y) 0<^<y (D23) 
vi{y)V2{x), y <x <ds ' 



v\{x) = coshfcg.x + 

ksKs 

V2{x^ = cosh fcs(a; — c?s). 



■ sinh ksx, (D24) 
(D25) 



Now, the solution of Eq. (D19) is 

F^^{x,iL0n)=2 dyG{x,y,iu;n)A{y). (D26) 
Jo 

Since the order parameter A{x) is determined by the 
symmetric part of the anomalous function F^ , the self- 
consistency equation comes from 

A{x) = wTcXg{ep) ^ F/(a;, iw„). (D27) 

Uln>0 

Substituting (D26) into (D27) gives the self- consistency 
equation, 

A{x)=2TrTcXg{ep) / dyG{x,y,iuJn)A{y). 

(D28) 

From this self-consistency equation, we can write down 
the equation for Tq with respect to the bulk transition 
temperature Tco- 



A{x)ln^ 



= 2nTc ^ / rfy 



Wn>0 



(D29) 

^ -G{x,y,i^n)\^{y) 



The above integral equation in (D29) can be reduced to 
an eigenvalue problem after we change the integration 
for the function G{x,y) over y into summation over dis- 
crete y. The integral domain y is divided into N grids, 
and we change the integration over y into the following 
summation to obtain Eq. (24) in Sec. 11. 



TV 



An In = KnmK. 



where we define 



(D30) 



(D31) 



a;„>0 



Knm = 2txTc V ^ - Gnn^iliOn) (D32) 

(D33) 



Gnmiioj) = G [ ^n,j^m,iuJri 



The Tc is determined by the smallest eigenvalue which 
gives the largest Tc- 
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